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ABSTRACT 

The  diffrectioa  of  plene  pulses  by  a  parabolic  cylinder 

pnd  by  e  paraboloid  of  revolution  is  investigated  for  both  the 

boundary  conditions  u  «=  0  and  *^  =  0.   In  all  esses  the  pulse  is 

on 

assumed  to  be  incident  on  the  vertex  of  the  diffracting  obstacle 
in  a  direction  pnrallel  to  its  axis.  The  problems  ere  treated  by 
first  sepereting  varipbles  in  the  time  dependent  wRve  equation  end 
then  solving  the  resulting  integral  eqiaations  by  meens  of  a  Laplace 
trnnsform.   In  eech  cese  a  Volterre  integral  equetion  for  the  solution 
is  obtained,  which  cen  be  solved  by  the  usviel  iteration  process.   It 
is  shown  that  when  the  pfirabolic  cylinder  degeneretes  into  a  half  plane, 
the  solution  reduces  to  that  obtained  for  the  corresponding  problem  by 
other  methods. 
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Ihtroductioa 

For  m?themstlcsl  es  well  83  physical  reesons  it  seems  convenient 
to  divide  dif fraction  problems  that  arise  in  the  theories  of  wpve 
propagation  governed  by  the  wave  eqioation 


^""bst 


c 

into  two  classes:   time  harmonic  problems  in  which  incident  tiiie 
harmonic  wrve  trains  are  considered^ rnd  problems  in  which  the  dependence 
upon  time  of  the  functions  involved  is  not  periodic.   In  the  latter  class 
a  typical  case  is  that  of  the  incident  pxilae  or  transient  wave.  A  pxilse 
wave  may  be  defined  for  the  purposes  of  this  article  as  a  moving  discon- 
tinuity surface  behind  which  there  is  a  prescribed  field  varying  in  space 
and  time  in  such  ?  way  as  to  satisfy  the  wave  equation,  and  in  front  of 
which  the  field  is  zero.  A  diffraction  problem  arises  when  sn  obstacle 
is  assumed  to  exist  in  the  path  of  the  moving  pulse  wave  and  the  field 
for  ell  time  is  constr?ined  to  satisfy  certain  homogeneous  boundary  con- 
ditions on  the  surface  of  the  obstacle.   In  addition,  the  physically 
reasonable  condition  is  usually  imposed,  at  least  by  implication,  that 
the  total  field  be  a  limit  of  continuous  solutions  of  the  wave  eq-uation 
and  that  a  first  time  derivative  of  the  field  be  a  limit  of  the  continuous 
first  time  derivatives  of  continuously  differentiable  solutions  of  the 
wave  equation.    The  diffraction  problem  then  is  to  determine  the  motion 
of  the  discontinuity  surface  and  the  field  generated  for  all  time  and  at 

ell  points  of  space  where  the  field  is  continuous. 

C7] 
In  teport  No.  SM-^3  the  problem  of  determining  the  field  due  to  an 

electromagnetic  pulse  wave  of  arbitrary  shape  incident  on  a  perfectly  con- 
ducting infinite  wedge  such  that  the  electric  vector  or  the  magnetic  vector 
is  polarized  parallel  to  the  generators  of  the  wedge  was  treated.   In  this 
report  we  shall  consider  the  diffraction  of  a  plane  pulse  wave  incident  on 
the  vertex  of  a  perfectly  conducting  parabolic  cylinder  in  a  direction  parallel 
to  the  axis  of  the  parabolic  cross  section  of  the  cylinder.  We  shall  again 


♦  For  a  discussion  of  the  implications  of  such  a  requirement  for  a  solution 
of  a  hyperbolic  partial  differential  equation,  the  conditions  which  the  dis- 
continuity of  the  solution  must  satisfy  as  well  as  the  manner  in  which  it 
must  propagate,  see  [l]  .  ?.  356  ff.  and  p.  470,  and  [k].     For  a  similar  dis- 
cussion of  a  discontfrmous  solution  of  Maxwell's  equations  in  electromagnetic 
theory,  see  [j^  and  l3[. 
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consider  two  cases:  one  In  which  the  electric  vector  of  the  incident 
field  is  essvimed  to  "be   polarized  in  s  direction  parallel  to  the  generators 
of  the  peraholic  cylinder,  end  ?nother  in  which  the  inegnetic  vector  of 
the  incident  field  is  esstuned  to  have  that  polarizetion.   In  these  cases 
the  prohlem  of  determining  the  field  reduces  to  solving  a  boundary  value 
problem  for  the  scalar  wave  equation 

XX   yy   g2    tt 

where,  for  case  1,  U  =  0  on  the  peraholic  cylinder  for  ell  time,  end  for 
case  2,  *—  =  0  on  the  paraholic  cylinder  for  all  time;  *—  refers  to  a 
derivative  taken  in  a  direction  normal  to  the  s\irfece  of  the  parabolic 
cylinder. 

We  shall  also  consider,  for  the  sake  of  completeness,  the  analogous 
boundary  value  problems  for  a  plane  pulse  similarly  incident  on  a  paraboloid 
of  revolution. 

Both  of  these  problems  in  the  case  of  sound  wave  propagation  have 
been  treated  by  Frledlander  '-^  and  by  Chester  ^  "^  .     Chester  obtained 
solutions  in  the  form  of  integrals  which  are  not  generally  suitable  for 
numerical  computation.  Frledlander  reduced  the  problem  to  one  for  a 
Volterre  integral  equation,  the  solution  of  which  can  be  obtained  by  an 
iteration  proced\irf'.   In  this  report  we  have  obtained  results  essentially 
equivalent  to  those  of  the  above  authors  by  a  different  method  which  is 
analogous  to  the  method  used  in  Report  No.  EM-^3  for  the  infinite  wedge. 
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X.  The  Parabolic  Cylinder 

1,1  Formulation  of  the  Characteristic  Boundary  Value  Problem 

We  shall  be  concerned  in  this  section  with  the  field  resulting 
from  the  diffraction  of  a  uniform  plane  pulse  wave  incident  on  the  con- 
vex side  of  a  parabolic  cylinder  and  traveling  at  the  propagation  speed 
1  in  e  direction  parallel  to  the  exis  of  the  perebolic  cross  section. 
We  shall  have  to  deal  with  the  wave  equation 

^xx*V-^tt"°'  (1-11^ 

The  incident  field  is   t^ken  to  "be  1  on  the  initial  side  of  the  plane 
discontinuity  surface  which  characterizes  the  incident  pulse  wave,  and  it 
is  taken  to  be   zero  on  the  other  side. 

The  equation  of  the  parabolic  cylinder  in  Cartesian  coordinates 
is  taken  to  be 

x^  -  Up(p-y)  (1,12) 

where  p  is  the  focal  length  of  the  parabolic  cross  sectiono 

In  (x,y,t)  space  the  incident  wave  front  is  a  plane  slanted  at 
U5  degrees  from  the  horizontal  (x,y)  plane.  Through  the  curve  of  inter- 
section of  the  incident  wave  front  and  the  parabolic  cylinder  in  (x,y,t) 
space  there  passes  another  characteristic  surface  of  (l,ll)«  This  second 
characteristic  surface,  representing  the  space-time  history  of  the  reflected 
wave  front,  is  a  characteristic  cone  whose  vertex  lies  on  the  focal  line  of 
the  cylinder.  The  fact  that  the  cone  passes  through  the  curve  of  intersec- 
tion of  the  incident  plane  and  the  parabolic  cylinder  can  be  seen  easily  from 
the  relations  t=p-y  for  the  incident  wave,  x  =^p(p-y)  for  the  parabolic 


/x2  2^ 


cylinder,  and  t+p=  Vx  +y  for  the  characteristic  conej  we  pssume  thet  the 
incident  wave  first  strikes  the  cylinder  at  time  t-O, 

In  place  of  Cartesian  coordinates  we  shall  use  a  system  such  thsL 
the  characteristic  cone  representing  the  reflected  wave  front  in  (x,y;t) 
space  is  a  constant  coordinate  surface  and  the  parabolic  cylinder 
representing  the  diffracting  obstacle  is  another.  If  t«0  is  the  time  at 
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VThich  the  incident  wave   strikes  the  vertex  of  the  parabola,  we  have  then 


t  -  I  U^*  vi^)  -((|+P)  (|  -  V^xV-(t+p) 

X  «  C^  "l^  -  Vx^+y^+y.  (1.13) 

C  and    T|      are  the  usual  parabolic  coordinates,    7^    takes  only 
positive  values,  but  4  i"ay  be  either  positive  or  negative.     In  these 
coordinates  the  cylinder  (1.12)   is  given  by 

y[    =V2p.  (i.ia) 

Employing  the  transformation  (l«13),  equation  (l.llj  becomes 

^C4  *  %  *2^^44  *   2^^<^^  *2^4  =  °*  (1.15) 

Now  if  the  reflected  wave  front  is  represented  by  the  expression 
T(x,y;=t,  then  T(x,y)  must  satisfy  the  eiconal  equation 

(VT)^-l. 

Since  the  characteristic  cone  whose  spatial  cross  sectione  are  the  various 

positions  of  the  reflected  wave  front  can  be  represented  by  t+p  =  yx  +  y  , 

this  requirement  has  been  met.  The  discontinuity  on  the  reflected  wave 

front  must  have  the  form^''''-^-l'  L^J'L^-1 

1  6 

[ul  -  D  e^o(^Y)ds, 

where  U  is  some  initial  value  on  each  ray.  If  we  use  polar  coordinates  (rj0) 
instead  of  Cartesian  coordinates,  we  can  write  for  the  discontinuity 

where  f(9)  is  a  function  to  be  determined  by  the  boundary  condition  on  the 
parabolic  cylinder. 

In  polar  coordinates  the  equation  (I.IU)  for  the  parabolic  cylin- 
der becomes 

^  2p(l-sinQ)   ^   2p 

2^       1+  sine  • 
cos  0 
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The  first  boundary  value  problem  is  that  in  which  u  »  0  on  the  dif- 
fracting obstacle,  or  u{K,V2p,   0)  =  0  .  Then  we  have   ^^^^  »  f(e)  \/J*^^  .  i. 

Thus  f(e)  =  /-,  ^  P  A   •  It  follows  from  (1.16)  and  the  fact  that  the 

incident  fielci  is  one  behind  the  incident  wpve  front  that  the  inhomogeneous 
boundrry  condition  on  the  inside  of  the  reflected  wave  front  is 


u(g,>;,o)»i--/  ,^^^P^^.  ; 


r(l+sin  0) 

9^  «  0  is  the  equation  of  the  characteristic  cone  which  represents  the 
reflected  wave  front.  In  parabolic  coordinates  this  condition  is 

u(C,n»0)  -  1  -  ^^  .  (1.17) 

The  first  characteristic  boundary  value  problem  can  now  be  stated: 

we  must  find  a  solution  u(4,>],9f)  of  (1.15)  such  that  u(4,V2,,0)  «  1  -V^^ 

and  u(C,/^  ,0)  "  0  in  the  region  determined  by  -oo  <  4  <oo,  V^ -'^^oo  , 
and  -co  <  (?  <  0. 

The  boxindary  condition  —  =  0  implies  that  the  field  inside  the  wave 

front  at  the  parabolic  cylinder  is  twice  the  incident  field  at  that  point,  or  2, 
Thus,  in  this  case  we  have 

u(4,'>l,0)  =  1  +  >^  .  (1.18) 

In  the  coordinate  system  given  by  (1.13)  we  have  on  the  parabolic 
cylinder 

Thus  the  second  characteristic  boundary  value  problem  is  the  following:  we  must 
find  a  solution  u(4,>z,9')  of  (1.15)  such  that  u(C,>2,0)  »  1  +i2E_  and 


^|i  .^^ii 


ri=/2?   ^^S./2i 


0 


in  the  region  determined  by  -co  <  C  <oo   ,    y^  ^V[<oo,   and  -oo  <  ^  ^  0. 
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1,2  Solution  of  the  First  Characteristic  Boundary  Value  Problem 
The  boundary  condition  (1.17)  involves  the  variable  t^  alone. 
The  boundary  condition  on  the  parabolic  cylinder  is  zero.  Thus,  it  seems 
reasonable  to  seek  a  solution  which  is  independent  of  4«  If  ve  assume  that 
our  solution  is  independent  of  4>  the  equation  (1.15)  takes  the  form 

We  must  find  u(i^,^)  satisfying  (1.21)  in  the  region  -co  <  ^  ^  0   ,/2p  ^>l<  oo 
with  u(v/^,^)  =  0  and  u(Ti,0)  -  1  - 
(1.21)  separates  into 


^-kf(i2)=0 

and 

(1.22) 

d\(n)      dH  (yj) 

d  V|  ^ 

where  k  is  the  sepsrstion  coastent,  J  is  a  function  of  ^  elone,  and  H. 
is  8  function  of  ^   elone.  We  pssume  s  solution  which  is  a  linear  com- 
bination of  solutions  of  the  equations  (1*22)  in  product  pairs: 

u()2,9>)  -/  A  (k)  Hj^(>2)e'^  dk  ,  (1.23) 

where  A(k)  and  the  contour  L  are  to  be  determined. 
Since  ui]/^,(jS)   »  0,  .t  is  required  that  \(V^)   =  0.   Moreover,  from 

u(i7,0)  ■=  1  -  ^  we  have 

/A(k)  H^(n)dk  -  1  -  ^  .  (1.2U) 

The   solution  of  the  second  equation  in  (1.22)  which  satisfies 

H.  (^)  =  0  is 

2     /-»!  2 

We  now  have,  combining  (1.25)  and  (1.23)  and  introducing  a  factor  2  for  future 
convenience , 

u(n,i2)  =  2/  A{k)e^^'^^^^  /"^   e-^^  dTdk.  (1.26) 
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We  require,  according  to  (1.2U), 

u(n,0)  -  2/    A(k)e^'iy_ 


2    r^    ^                       V^ 
^     '        e^^    dtdk  -  1 

2 

or, setting    >]     »  x,  - 


/  A(k)e^^    /    ^^dTdk  =  1  -V^  . 
L  2p    /f  ^ 


Then,  changing  variables,  >%  have 


x-2p         _^^  ^ 

/  A(k)e^(^-2PV  -^ dTdk  .  1  -/1e     , 

L  o  /r+^p  ^ 


and  setting  x-2p-  z,    we  find 

/A(k)e^V'7^     dTdk-l-/^     . 

This  suggests   that  we  might  make  L  the  contour   (c  -  ioD,  c  +  iOD  } 
•ased  in  the   inversion  theorem  for  the  Leplace   trsnsform.      Continuing 
formally,  we  heve 

z 


o     v  +2p       L  ^ 


and  we  observe  that  the  left  side  is  a  convolution  integral  for  the  Laplace 

transform.  Thus,   taking  Laplace  transforms  of  both  sides,  we  have 

00         -kT                                   00       -kT 
A(k)/     — ^ dT-i-v^/     -2 dT.  (1.27) 

0     /r+2p  0     /f+^ 


Finally, 


A(k) i -V^     . 

k/  -2— idr 

o    y/T>2p 
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00 


-kT 


We  introduce  the  notation  Lj^  {f  (T)^  «  f   f  (T)e         dT      for  the 

o 

Laplace  transform.     Then,   substituting  (1.2?)   in  (1.26)  with  the  appropriate 

change  of  variables,  we  obtein 


u(r|,5»)  =./  A(k)e*^(^-^y"?:^drdk-/J         ,\      , 
L  2p  v^  L      kU 


-^ 


^k(x^?)y.^  ll^Tdk 
2p/F 


-e-Lj-^ 


V^^^ 


dk 


./e;<^^-P),.^^^^K(x-,-.p)        ,^ 


/r^ 


dk  - 


^■"^f-^ 


kL 


(•r^j 


dk 


^L  ^l/pTTj 


dk   . 


The  first,  second,  and  fourth  integrals  are  simply  inversions  of  the 
Laplace  transforms  which  appear  in  their  integrands.  The  second  and  fourth 
integrals  cancel,  and  the  first  integral  has  the  value  1  •  Thus 

u(72,?)  -  1  -  K-?'),  (1.28) 


•vrtiere 


or 


e^L. 


i(y) 


1^1 


L      k 


W 


dk. 


(1.29) 


Vf+2p  J 
We  now  take  the  Laplace  transform  of  both  sides  of  (1.29)  end  obtain 


\ 


{....j-s 


■^1 


M 


/r-»-  X 


VT'SpJ 


^#iKH-^wM^} 
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Thus,  by  the  faltung  theorem , 


/    i(r)dr      _  /    dr 

o      l/y+2p-T  o    /r+x 

Differentiating  both  sides  with  respect  to  y  we  have 

0     (y+2p-T)  ^ 

Instead  of  attempting  to  justify  the  formal  derivation  of  (1.210),  one 
can  show  directly  that  u(>2,^),  given  by  (1,28),  with  l{-^)  obtained  by  solving 
the  integral  equation  (1.210),   satisfies  all  the  conditions  of  the  problem  of 
which  u(l2,9')  is  a  solution. 

If  we  set  y=0,  we  obtain  1(0)   =  y  ^  =  ^^     .   This  verifies  the  fact 

V2d 
that  u(r[,0)   »  1  -  ^L^  .     To  obtain  the  other  boundary  condition  u(-/2p,9f)  -  0 

vre  integrate    (1.210)   from  0  to  y  on  both  sides  of  the  equation  and  obtain 
again 

/      KT)  dr       ^   /      dr 
o       /y+2p  -T  o    /T+x 

which  cpn  be  written 

/  i(t)dr        /      dr 


Putting  "h  =  Y2p     we  obtain 

y  dY  -  0. 

o    yy  +2p-T 
Dlfferen tip ting  this   integral  with  respect  to  y  we  ostein 

/2r  o     (y+2p-r)^/'^ 

Hence   l(0)  =  1.     How,    integrating  by  parts,  we  obtain 

0      /y4.2p-T 
which  is  of  the  same  form  as   the  original   integral,    except  thet  we  have   I'(^) 
insteed  of  l('^)-l.      In  the  seme  manner  we  cen  show  that  all 
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derive tlTss  of  l(y)  ere  zero.  Since  we  shell  coneider  only  e  reguler 

solution  for  l(y),  it  follows  that  l(y)  =  1  when  7]^=  ^/^  J  i.e.,  u(  V2p,  (>)   =  0, 

end  our  homogeneous  boundary  condition  is  setisfied, 

Similerly,  since  the  right  sids  of  (1.210)  satisfies  (1.21)  and 
"becpuse  of  the  houndpry  condition  on  I  ft  zero,  we  cen  show 


/ 


T372 


dr=  0. 


o   (-^  +  2p-r)- 

From  the  regularity  of  I  it  follows  that  I  satisfies(l,21) . 
We  can  solve  (l,210)  by  iteration,  setting 


■■1       ''  x+y        '^  2    J 


io(r)dr 


"^  0       (y+2p-T0^/^ 


i-,(r)dr 


^^'^^^^/^^i;:^'-' 


I  =/I1p  ./I/  -^ 

n       'x+y        '    2  J        f 


i„  _L(r)dr 


-y     ''i     iy.2v-r)y^  • 


r-       y   i      it;-!    _(,t;        ,17 
ve  have  then  |  I^^.lJ  =  /|  |  /  ^^^^^  _^^,%'---  ^^^  1^  ^/   |I„     -I„.i  I  «• 


If  vge   assume  y^  0,   as  it  is  for  all  points  exterior  to  the  parabolic  cylinder. 


and  set  I    =0,  we  have 

0 


n 


Thus  I  =  y~  (I.  T-I.)  approaches  a  limit  dominated  by  — ^  5~  ^  =  TP  ® 

1=0  ^      n=0  (Up)        '■ 

which  converges  rapidly  as  long  as  y  is  not  >>p.  Moreover,  since  the  I  are  all 

regular  and  uniformly  bounded  in  any  closed  finite  region  for  y»0,the  limit  I  is 

a  regular  function  in  any  such  region. 


-li- 


lt can  be  seen  that  the  iterations  are  not  difficult  to  carry  out, 
since  the  integrations  involved  are  elementary.     For  example,  we  have 

T   (y)  ./IE  ,   and  I,(y)  -       ^^V^^^       .     ^  ^P         -fill  ^  yTIl 
■^1'"^'^       '   x+y   '  ■^2^-'''  2p+x+y  2p  +x+y    '    2p+y      'x+y     * 

so  that  for  two  iterations. 


121 /IST.    V2p(>;^-g 


u(r^,/)-l  +  ^ /-f^    -    y^P^%-^^       .         (1.211) 

(2p  +  »i  -9)/2i^  >2   -^  2p+>f-  91 

This   ia  8  good  epproximetion,   et   least  if  -^  <  4p,    i.e.,    if  the  field 
is  meesured  et  s  point  not  too  lon^  efter  t  w?ve  front  hfs  passed. 

For  the  cpse  p*=0,   the  parabolic  cylinder  degenerftes   into 
the  hrlf  plpne  Pnd   (1.29)  "becomes 

^/T     L    /F     o       ^/t+  X  "0      y^^Ir  y?T3i  "  y 

where  we  a?ain  make  use  of  the  faltung  theorem. 
This  result  agrees  with  that  obtained  >hen  the  general  solution  to  the  problem 

\i\ 

of  a  pulse  diffrated  by  a  conducting  wedge  treated  in  EH-U3'-  -^  is  specialized 
to  the  Case  considered  here.  As  expected  this  result  also  agrees  with  that 
obtained  by  Keller  and  Blank, for  their  case,  by  a  different  method,'--' 

den  we  have 

Z^  (1.213) 


u(^,<!f)  -n^^^ 


1,3  Solution  of  the  Second  Characteristic  Boundary  Value  Problem 
We  start   as  before  vdth  the  assumption 

^{\y9)  -/   A(k)e"^  Hj^()Z)   dk. 

Ju 

According  to  (1.19)  we  must  have 


^1  »y  A(k)e"^       —^ ky^Hj^(>Z)  dk  -  0.  according  to 


du 


1=)^ 


We  need,  therefore,   a  solution  ofs^the   second  equation  of  (1.22)  vhich 


satisfies 


dH.(Y?) 


V>?)] 


-knv^?)       _  "0  • 


^n     "^'*k-'j^.^ 
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The  solution  we  desire  is 


.2 


\(>2)-e^^     (1-/5?  k  e^PV'^    e-^^    dX ) ,  (1.31) 

and  we  assume 

L  y25 

which  SEtisfies   the  boujidnry  condition 

2  n  2 

u(rj,0)  -/  ACkje^"^     (1-/5?  ke^P^/       e"^'^    dfjdk. 
L  /2p 

We  observe   that  the  constant    1     satisfies  the  differential  equation  (1.21)   and 

the  boundary  condition  -s—    |  =*  0    . 

Thus  we  can  write 

u(n,9i)  -  l+/A(k)e^^'^^-<^^   (1-/2?  k  e^P'^/     e"^''    dr)dk  (1.32) 

L  /2p 

and 

.u(>2,0).l  =/A(k)e^^V-/5?  ke^P*^/"^  e'^^'^dTO  dk  =    ^  .  (1.33) 

L  /2p  ^ 

If  we  let  z  =^    -2p  pnd  B(k)  =  A(k)e^'^  ,   this  lest  equfstion  is  equivalent  to 


z         -\Ci^ 

r       e 

L  '   '^''l  "o  yr +2p 


/  B(k)e^'=  dk  -  y|/l^B(k)e^V'  ^^^  <i^  ^k  -  /^ 


If  we  proceed  formally  as  before  to  solve  for  B(k),  we  obtain 

/  B(.)e'<-di<  -  /f/'  -i^  /kB(.)  e-t-^dl-  -  ^^ 

L  o      /r  +  2p         L  ^ 

and  finally  00       .jjT 

/2?/    -2 dr 


However,    if  we  were   to   continue  rs  before   to   obtain  e  Volterre    integrpl 
equation  for  our  solution,  using  this  vplue  of  B(k),  we  would  see   th?t  our 
boundpry  condition   (1.33)  will  not  be  satisfied.      It  is  easy  to  see,  nevertheless, 
thpt  P  vrlue  for  J(k)  which  does  work  is   just  one-hrlf  thft   obteined  formelly. 
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The  incorrect  value  of  the  formal  result  is  evidently  due  to  the  fact  that 
the  interchange  of  integral  signs  in  our  procedure  is  not  legitimate  in  this 
case .  ^-'e  proceed  therefore  with  the  value 


00  -k 

e 


B(i<)   =         ^   -^o       /rT2f 


dr 


o    »^+2d 


and 


h 


u(^p)  =  i^/iy  — M 


-k/i 


h 


r-   r     -kT 


o  /r+2p 


By  vrritingy 


z  -kT 

e 


dT 


yr+2p    J 

CO  00 


-kT:' 


(1.3U) 


dT;6k^^-^)dk.  (1.35) 


d"^      ,  we  can  split  the 


o    yr+2p  ■'  o     "  z     /T+2p 

integral  on  the  right   side  of  (1,35)   into  two  parts,   one   of  which  can  be  evalu- 


ated iramedi 


iately:  /|   /  e^^^"^^   L    [— ^^  /dk  +  l(-^,z)   -  /|         \  +  l(-9r,z)    , 


where 


kL, 


^ 


i(y,z)  .  H    /   tei^pj /  )/rT2^^J    ^ky 


dk. 


(1.36) 


/l^+2p 


We  can  find  a  Volterra  integral  equation  for  l(y)  just  as  before.  First,  it  is 
necessary  to  integrate  by  parts  in  the  numerator  and  denominator  of  the  integrand 
of  (1,36),  so  that  the  expression  for  l(y)  becomes 


:(y)  -  /I. 


\ 


2/ 


y/r  -*•  z-t-2P 


1  *  /f  ^  [777-0/2 


e'^dk-E 


lie 


\i^H^^ 


}(T*2p)3/2 


^^'^'^f^^^J 


:^^^2pJ  gky^k. 
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We  now  take  a  Laplace  transform  of  both  sides  of  the  equation,  clear  fractions, 
and  take  the  inverse  Laplace  transform,  applying  the  faltung  theorem  as  before; 
we  ol)tPin 

^     o     (y-r+2p)^^^       /y+z+2p       o  (y  .r+2p)^/'^  yr*  z+2p 
thft  is, 

the  whole  solution. we  remember,  is  given  by 


u(yi,(2)  -  1  +  /   ^    +  I(-9I)  .  (1.38) 

We  observe  that  according  to   (1.37). 

1(0)  =     ^/|. 
so  that  along  with  (1.38)  we  have  u(>/,0)  =  1  +   -^^-^  ^  as  required. 

To  check  the  boundary  condition     Si  ■  0    we  require   the 

expresaiort     —  I  =  v  .  which  is   the  solution  of  the  integral  equation 

^'^    ^=/2? 
Obtained  from  (1.37)   oy  differentieting  end  setting    ^  = /^  .  nemely 

^■^o      (2p-^T)5/2        (^p_j^)2  2(^^)-'/'^ 

It  cen  be  verified  that  the  only  solution  of  this  equetion  is 

^  =  I  (2p-^)"^/^ 

If  we  compute   the  normel  derivrtive  of  u  from  equetion  (1.38)  for   ^  =  y^ 
pnd  insert  the  sbove  vpl->je  for  vf,  we  find    LL:.  t  r—  |  =  0,  as  required. 

By  p  procedure  similpr  to  the  one  employed  in  Section  1.2  it   is  e  mstter  of 
straightforward  calculations  to  show  that  u{rj  ,9^)  given  by  (1.38)   and  (1.37) 
satisfies  the  differential  equation  (1.21). 

We  can  solve   (l,37)   as  before  by  iteration.     The  result  converges 

rapidly  toward  the   correct  solution, provided  again  that     -l2_  is  not  too  large. 

*  Up 

For  exemple,    zero-order  epproiimetion  in  this  esse  will  be 

wvj,^)  .l./ZlZr.    E Pl^T       .    _^ .(1.39) 

2{ri~9)  nV^vT  >/(2p-9+r)         y>^2_^  (2p.0+x^2) 
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2.   The  Paraboloid  of  Revolution 

2.1  Formulation  of  the  Characteristic  Boundary  Value  Problem 

In  this  case  the  parabolic  cylinder  is  replaced  by  a  paraboloid  of 
revolution  given  by 

x^  +  y^  =  Up(p-z)  ,  (2.11) 


where  p  is  again  the  focal  length.       The  wave  equation,    similerly  in  three 
dimensions,   is  ^  ,  ^ 

In  place  of  the  Cartesian  coordinates  (x,y,z,t)  we  shall  again  use 
a  coordinate  system  adapted  to  the   special  shape  of  our  obstacle,  the  para- 
boloid of  revclution,  and  to  the  characteristic  surface  representing  the  re- 
flected wave  front.     Wa  write,  therefore. 


P  -  /x^+y^+z^    -  (t+p)  ,  t  =  i  iK^^yf)  -(g+p) 


(2.13) 


\        -  yx   +y  +z      +  z  ,  X  =  4 17  cos  Y 

2      n — 5 — ? 

?"V«+y+z-z,  y  =  4V)sin  T 

T    -  tan"^^  ,  2  =  \{^-h  . 

In  these  coordinates  the  paraboloid  of  revolution  becomes 

>?  -  /2p   .  (2.11i) 

Becpuse  of  symmetry  we  may  assume  the  wave  equation   (2.12)   to  be   independent  of  y^\ 
then 

We  can  show,    In  a  manner  completely  analogouB   to  that  employed  in 
the  case  of   the  parabolic  cylinder,    that  the  reflected  wave  front,  which  is  at 
each  instant  of  time  a  cross  section  of   the  four-dinensional  characteristic 
hyp«Tcone  whose  vertex  is  at  the  focus   of  the  paraboloid  of  revolution,    is  a 
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sphere  whose  center  is  at   the  focus.     This  time  the  discontinuity  conditions 
show  that   the  discontinuity  across  the   spherical  wave  front  reflected  from  the 
paraboloid  varies  in  such  a  manner  along  a  ray  that  it  has  the   form 

fuj   .     f(Ml  ,  (2.16) 

where  {r,Q, (jl)   are  ordinary  polar  coordinates  with  the  coordinate  system  center 
at  the  focus.   The  function  f(r,0)  can  be  determined  'by   inserting  in  (2,16) 
the  velues  of  ^ii]  pnd  r  on  the  pereboloid.  The  equation  of  the  parabloid,  in 
polar  coordinates,  Is  given  by 

Fnd  for  the  boundf ry  condition  u  =  0,  [u]  =  1  on  the  surface.  It  follows 

1  +  cos  0      *- -'   r(l  +  cos  Qi).        In  pereboloidel  coordinptes 
[u]  =  =^  end  it  ag-tin  eppeers  that  [u]  depends  on  ^elone.  We  now  have 


u(>i,4,0)    =  1  -     ll     ,  (2.18) 


-^  +  Uu„  -  0    .  (2.19) 


and  we  can   assume  once  more  that  u  is  independent  of  C«     With  this  assumption 
(2.l5)  reduces  oo 

^vi>v    *   2v|u^^    ^  ^u^ 

The  first  boundary  value  problem  for  this  case  is  then  to  obtain  a 
solution  u(  >|,9)   of  (  2JL9)   in  the  region  defined  by  -oo  <  $2  <  0,  0  ^  T^  <oo 

such  that  u(>T  ,0)   =  1  -  ^  and  u(/2p,^)   -  0. 

The   second  boundary  value  problem  can  be   set  up  in  a  way  exactly 
analogous  to   the  second  boundary  value  problem  associated  with  the  parabolic 
cylinder.     The  boundary  condition  on  the  paraboloid  of  revolution  is 

3u    I  Two"     du  ^   3u  c\ 

By  the   same  argument  used  before  we   shall  have  1  +  — §  for  u  on  the   reflected 
wave  front.     We  must  determine  a  solution  u(>2  ,0)  of  (2.19)   in  the  region 

determined  by-oo<9f<0,0<V|  <oo   such  that     |^    |  =  0  and  u(>^,0)=l+  -|. 


is 


\{yi)   '  ^^^    J  V-     dT     .  (2.22) 
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2.2  Solution  of  the  First  Characteristic  Boundary  Value  Problem 
Equation  (2.19)    separates  into 

J(gf)   -  -k^(9J)  (2.21) 

^^  h|^(>?)   -  2k>^H^(>/)  +  i  Hj^(>?)   -  iac  Hj^(>/)   -  0    . 

The  solution  of  the  second  equation  of  (2.21)  vhich  vanishes  at  ^  =  /2p 

We  therefore  assume   thet 

u(jl,g)  =  2/ A(k)e^^'^  -^V        V-    dTdk  (2.23) 

pnd,    in  ©ccordrnce  with  the  boundary  condition  on  the   reflected  wgve  front, 

that 

/2    /"^     -kf 
A(k)e'''^    /        ^-r-  dYdk  -  1  -  ^    .     (2.2U) 

The  integral  equation  (2.2^+)  cfin  he  solved  in  exactly  the  same  way  as  the  corres- 
ponding integral  equation  associated  with  the  parabolic  cylinder.  Writing 

2 

z  «>2  -2p,  we  have 

/'A(k)e^^  /  -^^  dYdk-1--^  . 
J  ^  J        ^  +2p  z+2p 

L  o 

Continuing  formally,  we  find 

f     .=:^   f  A(k)e^(-^  dr  -  1  -  -i^ 
-/^  T  +2p  v/j^  '  '  z+2p 

from  which,  using  the  faitung  theorem,  we  obtain 

°°         -k^ 


1  -  2pk 


A(k)  -    „  °      \v^P .  (2.25) 


00  -k't' 

'-       e_ 


"/    -h^  '^ 
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On  substituting  this  result  into  (2.23)  we  obtain 


^klYTTp^     •  '' 


21^   -1^ 

®    dTdk 


°°   °°  -kT 


^dT^dk 


^<^'  -^-^^'a.  V:,[^5}««^^^-^>- 


^  --L 


or,  finally, 


u(y^,9r)  =  l-I(-i2),  (2.26) 

where     l(y)  -  /   ■■■  ^,   dk .  (2.2?) 

¥e  can  now  obtain  a  Volterra  integral  equation  for  l(y)  by  the  pro- 
cedure we  have  used  before.  Taking  the  Laplace  transform  of  both  sides  of 
(2.2?)  and  clearing  fractions,  we  obtain 

By  the  faltung  theorem,  then, 

y  (y.2p-r)  ^^   /  r-^r^2   ^ 
and  differentiating  with  respect  to  y  we  heve 

I(y).jp/  m^I-^       .       -a_   ,  (2.28) 
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From  (U.28)  we  have  l(o)  »  -S  which,  combined  with  (2,26),  shows 


^(yifO)   "1 ?  ,  ^^  desired. 


As  before  we  can  show  that  u(?7,^)    as  given  by  (2.26)   with  l(y)  de- 
fined by  (2.28)   satisfies  our  differential  equation  and  the  boiindary  condition 

u(y^,0)  =  0. 

We  nipy  obtain  I  by  iterntlng   (2.28)  fs  before.     The  solution  will 
pgrin  be  vfilid   in  the  neighborhood  of   the  reflected  weve  front,   with  the 
region  of  vrlidity  vprying   inversely  es   the   focpl   length  p. 


2,3  Solution  of  the  Second  Characteristic  Boundary  Value  Problem 

We  now  consider  the  boundary  condition     —  =  0  on  the  paraboloid 


of  revolution.     Once  again  we  have  ^ 


\-^ 


[^l*i|] 


)2«y^ 


0  . 


In  place  of  u(>^,0)   -  1 ^  we  have  u(i^,0)   =  1  +  -^    on  the  wave  front. 

The  solution  of  the   second  equation  of  (2.21)   required  in  this 


ce-se  is 


HJ>?) 


r 


n    -k^^ 


l-2pke^PV  V 


dr 


M 


Thus, 


u(>?,g)   =  1  W     A(k)e''^'^  -^^     l-2pk6^P^  / 


2pk   /  e  ^  ' 


T 


d^ 


dk         (2,31) 


pnd 


u(>^,0)   -  1  +y  A(k)   e*"^ 


l-2pke 


2pk 


■n   -kr2 


y^ 


r 


<rt 


-  =  % 


+  1 


thnt  is, 


/"A(k)   e^^^   L2pke2pl^/ 


^    -k 

e 


-k^ 


'T 


drt 


dk=  ^, 

1' 


(2,32) 


We  must  solve  for  A(k)  in  (2.32).  We  proceed  formplly  ps  in  the  previous  cpsea, 
setting  >]^-2p  »  z  and  B(k)  -  A(k)e  ^  and  changing  variables  in  the  integral 
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so   thpt 

/  B(k)e>=^  *-p/    ^J  ^(.)e«--'  ^  -    ^  . 
L  o  L 

whereupon  we  heve   the  forme 1  result 

00        -kr 

o 
But  in  this  cese,  es  in  a  previous  one,   this  vslue  of  B(k)   is   incorrect,  end 
for  the  spme  re?son.     However,   one-half  the  above  result  works,   as  we  can 
demons trrte:     We   set 


-    1   + 


iT^  *  P«  A 


thnt  is, 


u(>7,9f)  =  1  +    —^    +  pl(-9f)  ,  (2.33) 


wherf 


1     .1^ 


dk 


dk 


Then  takine  Laplace   transforms  and  applying  the   faltung  theorem 


we  obtain 

y 


"^o    (y-r+2pr       ^+y  o     (r+2p)^  (y  +  >?'^--l)     ' 

thet   is 


> 


,(y,  .jp/^(l)dl  -A. 


o   (y+2p-t)^       y+>^^  (2p+y)(>|^+y+2p) 


2piogj(^p;W^(/-y)]  (2.35) 


2p 
(>Ay+2p)^ 
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It  follows  that  l(o)  =  — ^     so  that  by  (U.33)  u(V7,0)  »  1  +  -^  as  required, 
¥e  can  show  thet 


k  t*i] 


^  I       -  |/2p  ^  +  ^  I       -  0  and  that  u 


and  that  u  satisfies  (2.19)  when  u  is  defined  by  (2.33)  snd  (2,35)  in  the 
seme  wsy  as  in  the  case  of  the  perebolic  cylinder.  The  equation  (2,35) 
cen  be  solved  by  iteration  es  before. 
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